Ferromagnetic Josephson junction with precessing magnetization 
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The Josephson current in a diffusive superconductor/ferromagnet/superconductor junction with 
precessing magnetization is calculated within the quasiclassical theory of superconductivity. When 
the junction is phase-biased, a stationary current (without a.c. component) can flow through it 
despite the non-equilibrium condition. A large critical current is predicted due to a dynamically 
induced long range triplet proximity effect. Such effect could be observed in a conventional hybrid 
device close to the ferromagnetic resonance. 

PACS numbers: 74.45.-l-c, 76.50.-|-g. 



The proximity effect in a ferromagnetic (F) metal in 
contact with a conventional superconductor (S) is usu- 
ally short ranged [l| . Indeed, the conversion from Cooper 
pairs in the superconductor to Andreev pairs in the ferro- 
magnet involves two electrons with opposite spins. They 
get quickly dephased due the large value of the exchange 
field acting on their spins. Yet, it was predicted that 
a long range proximity effect can take place when the 
magnetization rotates spatially [2] or when spin-flip pro- 
cesses take place [Sl in vicinity of the F/S interface. 
Indeed, Andreev pairs involving electrons with parallel 
spins are then created and can propagate on a long dis- 
tance. The long range proximity effect was invoked to 
explain the large Josephson current measured through 
half-metallic (HM) chromium oxide [4j as well as super- 
conducting phase-periodic conductance oscillations in an 
Andreev interferometer made of a long metallic wire with 
helimagnetic order [5|. However, the crossover from short 
range to long range proximity effect in the same device 
remains to be observed. Recently, it was suggested to 
measure the strong variations of the Josephson current 
through a trilayer ferromagnetic junction by tuning the 
relative orientations of the magnetizations [6| . 

In this article, we propose that the long range triplet 
proximity effect can be stimulated by varying in time 
(rather than in space) the orientation of the magnetiza- 
tion in the ferromagnet. Specifically, we consider a diffu- 
sive metallic S/F/S junction with dynamical precession 
of the magnetization, as shown in fig. [T] In the ferro- 
magnet, the conduction electrons feel the time-varying 
exchange field: 
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iF{t) = h F {sin 9 COS nt, Sin 6 sill nt, COS 9) (1) 

which is proportional to the magnetization and assumed 
to be spatially uniform. Here, r2/27r is the precession 
frequency around z-axis, is a tilt angle, and the ampli- 
tude of exchange field, hp, is constant at temperatures 
well below the Curie temperature. In the rotating frame, 
such precession can be viewed as a difference between 
spin-resolved chemical potentials in the leads, /is = fi^z. 
At finite tilt angle, hp has a component transverse to 
fis- Thus, it creates a non-equilibrium situation for the 



FIG. 1: Geometry of the ferromagnetic Josephson junction 
with precessing magnetization and configuration of d.c. and 
r.f. magnetic fields inducing the ferromagnetic resonance. 



conduction electrons. In particular, the device could act 
as a spin pump when the leads are in the normal state [7| . 
Moreover, in the superconducting state, the non colinear 
orientation of hp and fig generates a long range proxim- 
ity effect similar to the case studied in Ref. [6[. Despite 
the non-equilibrium condition, we find that a stationary 
(d.c. only) Josephson current flows through the junction 
when a phase difference is applied to the leads. In a long 
junction, its amplitude is comparable to that of a normal 
(N) metallic Josephson junction with the same length, 
provided that hfl is comparable with ksTc {Tc is the su- 
perconducting critical temperature of the leads) and 9 is 
large enough, while it would be exponentially suppressed 
in the absence of precession. 

Recently, a Josephson effect was also predicted in 
S/HM/S junctions when magnetic impurities located at 
the interfaces dynamically precess and mediate spin-flip 
tunneling processes for the conduction electrons [8|]. It 
was also suggested that this effect could be observed at 
the ferromagnetic resonance (FMR) . In the present work, 
we address the opposite case when the Fermi energy much 
exceeds the exchange field and the spin polarization in 
the ferromagnet is small. When the junction's length, 
L, much exceeds the ferromagnetic coherence length, 
^i? = y''hD/hp, where D is the diffusion constant, we 
will show that the Josephson current has two dominant 
contributions. The first one comes from long range triplet 



proximity effect. The second one arises from the inter- 
ference between short range and long range proximity 
effects out of equihbrium and it displays a non-analytical 
temperature dependence close to T^. Both have an oscil- 
latory behaviour, depending on the ratios between hQ, 
ksTc, and the Thouless energy Et = hD/L^. This pro- 
vides a mechanism for 7r-coupling [l| in long ferromag- 
netic Josephson junctions. In the following, we derive 
both components, discuss their properties and the con- 
ditions to observe them in conventional F/S structures. 
Within the quasiclassical theory of superconductivity, 
the current flowing through the junction is |9l.lld|: 
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Tr[r,(.9oV.9)^^(i,t)] 



(2) 



where the quasiclassical Green's function g{x,t,t') is a 
matrix in spin, Nambu, and Keldysh spaces and it obeys 
the Usadcl equation along the ferromagnetic layer: 

^iDVigo\/g) + [{idt + hF{t).(T)T,S{t-h)°g] =0. (3) 

Here, G is the conductance of the ferromagnetic layer, 
the spatial derivative is taken along x-axis, e is the el- 
ementary charge, at and Tj {i,j = x,y,z) are the Pauli 
matrices in spin and Nambu spaces, respectively, and o 
denotes the time convolution. (Units with h = ks ~ i 
are adopted from now.) Note also that the orbital ef- 
fect generated by the magnetization has been neglected 
in Eq. ([3]) as it is usually done [1[. Moreover, g obeys 
the properties: Tr g = and g o g = 1, and it has a 
triangular structure in Keldysh space with the retarded, 
advanced and Keldysh components: g^, g , and g^ , re- 
spectively. The normalisation condition is fulfilled pro- 
vided g^^"^ o ^^M = 1 and g^ = g^ ° f — f ° g"^ , where / 
is a distribution function matrix which obeys [/, r^] = 

In the present study, we assume that there is a good 
electric contact at the F/S interfaces and we neglect the 
inverse proximity effect in the leads. This yields the 
boundary conditions: 

gix = ±^,t,t')=g,^^^^2{t-t'). (4) 

Here, x is the phase difference between the leads. The 
quasiclassical Green's function in a superconductor with 
phase (p is defined in energy space by g^ L (e) = [—ieTz + 
Ar^/v/A2-(£±ir)2] and gf^^ = {gf:^-gt^)fT. where 
T0 — COS (j>Tx ~ sin (f>Ty, /T(e) = tanh(e/2T), and A is 
the conventional superconducting gap at temperature T. 
We also introduced a small, phenomenological depairing 
parameter F which may account for the current-induced 
depairing, as well as scattering on magnetic impurities in 
the leads. 

In order to determine the current through the junction, 
one has to solve Eqs. ©-(jj]). Despite a time- varying ex- 
change field in Usadel equation ([3]) , the problem is equiv- 



alent to a stationary (though non-equilibrium) one. In- 
deed, one may perform the unitary transformation: 



'^it,t')^VU{t)g{t,t')U{tyV^, 



(5) 



where U{t) = exp{iiltaz/2) transforms from the lab- 
oratory frame into a rotating frame and absorbs the 
time-dependent terms in Usadel equation ([3]), while V = 
exp(zcKTj,/2), with tana = hp sinO/ {hp cos 9 + fl/2), ro- 
tates the spin quantization axis. Specifically, V aligns the 
effective exchange field with amplitude J — [h'j^ sin^ 9 + 
{hp cos 9 + 17/2]^/^ along z-axis and rotates the preces- 
sion axis away from it. Eq. Q now takes a simple form 
in energy space: 

- iDV{g{x, e)Vg{x, e)) + [(e + Ja,)T,,g{x, e)] = 0. (6) 

The boundary conditions ([4]) yield: g{x = ±-§,£) = 
^5s>Tx/2(^+^'''2/2)^^ and the current ^ takes the form 



/=-§^ I deTT[Tz{gyg)^^ 



(7) 



Let us emphasize that the special time-dependence in 
Eq. II]) is crucial for the non-equilibrium problem ([3]) to 
be formulated as a stationary one, see Eq. ([6]). For a dif- 
ferent time-dependence and/or for applied bias voltage. 



one could use the procedure formulated in Ref [ll| to 
study numerically the conductance in S/N/S junctions. 

In spite of the above simplification, we are still dealing 
with a complicated non-linear equation ^. Now, we 
assume that the temperature is close to Tc, so that A ex 
\Tc{Tc — r)]^/^ is vanishingly small. Thus, we can look 
for a solution g = go+.9i +32 + - ■ • in the series expansion 
around the normal state solution, when A/Tc <^ 1. In 
zeroth order, one finds g^ = ±rz and 



/o = ./+-,/ 
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cos aa2 



(8) 

w here /± (e) = [frie + 0/2) ± frie - fJ/2)]/2 and q = 
\f2iJjD. We note that /o is diagonal in Nambu space 
and the non diagonal components in spin space only ap- 
pear at finite fi. In the first order in A, Eq. ([6]) yields: 
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which is solved with g^ = g^Q + Qi-cr, where: 
sinhfc(L/2 — x) 



(9) 



gix = -F^sina 



sinh kL 

smhk{L/2 + x) 



-^x/2 



.9fo ± giz 



sinh kL 

-- (i^f ± F^ cos a) 

smhp±{L/2 + x) 
sinh p±L 



^-x/2 



(10a) 



sinhp±(i/2 — x) 



sinhp-i-L 



-x/2 



-■^x/2 



(10b) 



and gfy 



R\U 



Here, k 



while .gf ^ -T^{9i) 
p± = ^-2i{e±J)/D, and F^ie) 



\F^{£ + 0/2) ± F '^js- ^/2)]/2, where F^(e) = 
i/\/ \J{e + ir^) — A^ [121 . We note that the component 
gix is long ranged, while gio and giz are short ranged. It 
is straightforward to check that the Keldysh component 
is solved with /i = and that the current ^ vanishes 
up to the first order in A. 

In the second order in A, the current is: 
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T, (2V/2 - g'^Vhgt 



-5R{[3fV5f - V5fgf]/o - [.9f]'V/o})] . (11 



(The terms with ^2 



of the identities: {Tz^g'-^'^} ± [di'^Y ~ coming from 



have been eliminated with help 

1^'^} ± (ff^ )^ 
the normalization condition g^ = 1.) The function /2 
vanishes at the boundaries with the leads and solves the 
differential equation: 

- 2iDV^f2 + 2[Jazj2] = -^D (Vpf gf V/o 

+V/o.g^V.gf + V(5f V/ogf )) • (12) 

Close to Tc, the Josephson relation remains sinusoidal, 
/ = IcSinx, and the critical current Ic can be obtained 
by evaluating Eq. (fTTjl . 

Before proceeding, let us make simplifying assump- 
tions. In usual ferromagnets, hp would exceed ft and 
Tc by several orders of magnitude. Therefore, a 
J 2:^ hp and p_ ~ q. We also assume that the junction is 
long {L 3> £,f)- Thus, we discard the exponentially small 
terms contributing to Ic and we find that it has two main 
terms: Ic = I^ + If . The first one comes from the long 
range triplet proximity effect (component glx)^ only: 
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(13) 



where uj ~ 
K^ = \/2{lu 



{2n + 1)ttT are Matsubara frequencies and 



iQ/2)/D. In particular, we get 
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hWTc 
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if < Tc < Et, 
if Tc < O < Et. 



Here, Ico = ttGA'^ /ieTc is the critical current for a S/N/S 
junction with Tc <^ Et at temperatures close to Tc- We 
notice that /* vanishes at O = ttTc- At larger frequency, 
it changes its sign and a 7r-coupling between the super- 
conducting leads is realized. When Q ~ Et, we also find 
an oscillatory behaviour of the frequency dependence of 
J* on the scale of Et (see fig. [U with an alternation of 
0- and TT-couplings. 




FIG. 2: Precession frequency dependence of long range (left) 
and "anomalous" (right) contributions to the critical current 
for ferromagnetic Josephson junctions with different lengths. 



In contrast to ferromagnetic Josephson junctions with 
spatial variation of the magnetization and S/HM/S junc- 
tions, there arises a mixing of short range and long 
range proximity effects through the non-equilibrium spin- 
dependent terms of the distribution function matrix in 
the system studied here. Most of the terms in Eq. pd|) 
reflecting such mixing are suppressed by the small factor 
£,F / L, when compared to /*. However, the contribution 
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G£,F sin cos 



def-\F2V^ 



kL 



sinh kL 



(15) 



coming from the "anomalous term" oc 5i*V/ogi in 
Eqs. (|lip - (fT^ needs more care. Indeed, for this term, a 
finite depairing parameter F is necessary to regularize the 
otherwise diverging integral over the energies. Assuming 
r <C A(T) and Vt <g; Et, we get in the logarithmic ap- 
proximation: 



GA 
32^ 



A n 

sin cos 6* In — tanh — . (16) 

t -i c 



This term vanishes hke {Tc - Tf^lTLiTc - T) close to 
Tc ■ Such non analytic temperature dependence was al- 
ready met in the field of transport in N/S junctions [13| . 
When il ~ Et, we also get from Eq. (fT51) an oscillatory 
frequency dependence of /^* (see fig. [5]). 

Let us now compare both main contributions p3p and 
P^ to Ic- They have quite different dependences on 
the device parameters. At low precession frequency com- 
(14) pared to Tc, Ic* scales linearly with O and dominates over 
/* which scales quadratically. On the other hand, at large 
frequency, /* takes over if A/Tc > (^f/^) ln(A/r). At 
low temperatures, we speculate that the above expres- 
sions still hold qualitatively. As the ratio A/Tc becomes 
of the order of unity, the critical current sd L ^ ^p is 
dominated by the long range triplet component only. 

Simplifying assumptions have been made to get analyt- 
ical expressions for the current. In principle, more realis- 
tic conditions could also be studied. The Usadcl equation 



^ may be solved numerically to get the current at all 
temperatures. Spin-flip diffusion in the ferromagnet and 
tunnel barriers at the F/S interfaces could easily be in- 
corporated in the theory. On the other hand, a finite 
spin polarization in the ferromagnet, as well as a mean 
free path comparable with the coherence length would re- 
quire to go beyond the quasiclassical theory for diffusive 
metals. Nevertheless, we expect that Eqs. (fT3 l) - (fT6| hold 
qualitatively beyond their strict range of validity. Thus, 
our findings are relevant to the study of F/S junctions 
with conventional metals. The quasiclassical theory has 
proved its usefulness for the prediction and quantitative 
analysis of several experiments in these systems [l| . 

The dynamically induced long range proximity effect 
studied in this work is expected to have special proper- 
ties. In particular, it would be of interest to character- 
ize the current-voltage characteristics in S/F and S/F/S 
junctions. We may expect an a.c. current response at fre- 
quencies mixing Q and the Joscphson frequency. Wc note 
this has been studied in the context of tunneling trans- 
port through a single magnetic impurity placed between 
two superconductors [1J| as well as S/HM/S junctions 
with extended interfaces [l5|. 

The long range Josephson current may be observed by 
erforming FMR experiments in the microwave regime 
in a planar S/F/S junction. This geometry was used 
in N/F/N junctions [171, ll8| to detect electrically the spin 
pumping [71 due to a precessing magnetization. A large 
resonance frequency, such that HQq ~ ksTc (IK corre- 
sponds to 30GHz), can be reached even by applying a 
moderate d.c. magnetic field Hq along the plane of the 
ferromagnet [16|]. A small transverse r.f. field. Hit, is 
used to induce the magnetization precession. The tilt 
angle 9 strongly depends on the precession frequency. At 
resonance, 9 can be estimated from the Landau-Lifshitz 
equations: 6 ^ fiHrf /ahflo, where /i is the Bohr magne- 
ton and a is the Gilbert damping parameter. In soft fcrro- 
magnets, a quite large 6 ~ 15° could be obtained. Away 
from resonance, 9 is much reduced. This would induce 
a strong frequency dependence, via angle 9, in the ex- 
pressions for the critical current derived above. We note 
that FMR was recently performed in an Nb/Permalloy 
bilayer [19[ . Below the superconducting critical tempera- 
ture, the reduction of the resonance width was observed. 
It was attributed to an efficient proximity effect leading 
to the reduced efficiency of spin-flip processes below Tc 
[20|. In particular, this shows that the proximity effect 
still persists at FMR. If the same conditions can be met 
in an S/F/S junction, we expect that the large critical 
current predicted in this work will be well measurable. 

In conclusion, wc have proposed a model of ferromag- 
netic Josephson junction with processing magnetization 
where a large current flows thanks to a dynamically in- 
duced long range proximity effect. The dependence of 
the current with the precession frequency shows an os- 



cillatory behavior. We have discussed the conditions for 
this effect to be observed under the condition of ferro- 
magnetic resonance in superconductor/ferromagnet de- 
vices with conventional metals. 
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